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A STABILITY STUDY AND MODAL ANALYSIS OF A 


SPACE -STATION— CENTRIFUGE CONFIGURATION 

By Robert K. Sleeper and Harold C. Lester 
Langley Research Center 

SUMMARY 

The planar motion of a centrifuge — space-station configuration has been examined 
to determine the types of instabilities that can occur. Although the study is devoted pri- 
marily to undamped symmetric configurations, the effects on stability of linear viscous 
damping and mass asymmetry were also examined. 

A linear analysis which was effected by assuming constant speed rotation showed 
that with little or no damping, unstable motion can occur in rotational speed ranges or 
regions representing two types of instabilities: (1) The first region was bounded by upper 
and lower critical speeds and in this region the motion consisted of a self-excited pure 
divergence similar to that found for rotating shafts with elliptical cross sections. All 
configurations exhibited this type of instability which was attributed to the bidirectional 
dynamic properties of the centrifuge. (2) Special attention has been directed toward a 
second type of instability which may be overlooked in an analysis but which can occur at 
rotational speeds above the critical speed; it consists of a self-excited oscillatory diver- 
gence similar to that experienced by helicopter pylon-rotor configurations. Configura- 
tions with a heavy centrifuge hub or light gondolas were more likely to have this type of 
instability. 

For undamped models a superposition of the complex eigenvectors is employed to 
represent the various types of real motion in the rotating reference frame and a special 
stability technique applicable to even-power polynomials is used to determine stability 
boundaries. 


INTRODUCTION 

A space station with an onboard rotating mechanism such as a centrifuge may 
exhibit instabilities such as those that have occurred with pylon-rotor configurations on 
helicopters in contact with the ground. For such configurations two types of instabilities 
have been found. 




One type of instability is similar to that exhibited by a rotating shaft with an ellipti- 
cal cross section, the critical speeds of which occur in pairs corresponding to the two 
principal stiffnesses. (See ref. 1.) Within the speed range bounded by the two critical 
speeds, the motion diverges in a nonoscillatory manner. The instability is classified as 
"self-excited" by virtue of the fact that the destabilizing force or forces are functions of 
the displacement vector and its derivatives. 

The helicopter analysis (refs. 2 and 3) also showed that a second type of instability 
characterized by an oscillatory divergent motion can occur at speeds above the critical - 
speed range. This instability is referred to as a "self-excited mechanical instability” in 
references 2 and 3; the term "flutter" which more usually denotes an aeroelastic phenom- 
enon (ref. 4) has been used by others (ref. 5). With either of these instabilities, kinetic 
energy of rotation is transformed into elastic energy. The energy transfer may be very 
rapid and may result in a catastrophic increase in the amplitude of the associated elastic 
deformation. 

In the present study, which represents an extension of a previous study (ref. 6), the 
stability characteristics of a space station with an onboard centrifuge rotating at constant 
speed are studied. The mathematical model includes linear viscous damping, but the 
effects of gravity have been neglected. Linear differential equations describing small 
motions in the plane of the centrifuge are derived and presented in a nondimensional form. 

Linear stability techniques are applied to determine and classify the regions of 
self-excited instability as a function of angular rotational speed for various combinations 
of the nondimensional system parameters. Furthermore, an eigenvalue -eigenvector 
solution is formulated by using an equivalent matrix of first-order equations of motion. 
Because of the lack of matrix symmetry, the eigenvectors are, in general, complex and 
hence classical vibration modes in the normal sense are not obtained. (See refs. 7 and 8.) 
However, it is demonstrated that the motion of the system can be described in terms of 
the complete set of eigenvectors for an arbitrary initial state vector. A Meerov-Fuller 
stability criterion (ref. 9) has been adapted to the analysis of the undamped system; it 
enabled a rapid determination of the parameter value combinations yielding unstable con- 
figurations. Although the effects of damping and mass asymmetries are demonstrated, 
most of the study is devoted to the undamped-symmetric model. 

SYMBOLS 

A-^n*)^^*),! coefficients of sixth-order stability polynomial 
Ag(fi*) J ( see eq * 

a^,a 2 axial displacements of gondolas due to lateral deflections (see fig. 2) 
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coefficient of viscous damping for hub support 
coefficient of viscous damping for centrifuge arms 
first-order dynamic matrix (see eq. (5)) 

uncoupled first-order dynamic matrix for undamped symmetric configura- 
tion (see eq. (14)) 

identity matrix 
hub support stiffness 

equivalent spring stiffness of centrifuge arms 
lengths of centrifuge arms 
total mass of centrifuge, mj + m 2 + 
mass of hub 

mass of space station (without centrifuge) 
gondola masses 

nondimensional Laplace variable, s/w e = a* + jo>* 
characteristic roots of system 

generalized coordinates associated with lateral modal displacements 
0j(p) and 0 2 (p)> respectively 

Laplace variable, a + 

time 

reference time 
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coordinate axes with origin at center of mass of composite system and which 
rotate with centrifuge (see fig. 2) 

space-fixed (inertial) axes with origin at center of mass of composite 
system (see fig. 2) 

state vector of four degree-of -freedom representation 

transformed components of state vector of four degree-of- 
freedom representation (see eqs. (2)) 


),x 3 ( ), 

),x 6 ( )^j 


displacement coordinates of centrifuge hub relative to space-station 
(less centrifuge) center of mass (see fig. 2) 

) transformed relative displacement coordinates of centrifuge hub 

displacement coordinates of space station (less centrifuge) from center 
of mass of composite system (see fig. 2) 

general complex response, Y 0 eP r 

complex amplitude 

redefined transformed state vector of four degree-of -freedom representation 
(see eq. (12)) 

),z 3 ( ) 

),h ( ) 
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components of transformed state vector (see eqs. (11)) 


distance from hub of effective hinge of centrifuge arm normalized to the 
arm length (see eq. (C6)) 

distance from hub of effective hinge of centrifuge arms normalized to arm 
lengths, Ll and L 2 , respectively 


rotation coordinate of centrifuge relative to space station (see fig. 2) 
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axial modal displacement factor of gondola masses due to lateral deflection 
viscous damping ratio of hub support relative to space station 
viscous damping ratio for centrifuge arms relative to hub 
diagonal eigenvalue matrix (see eq. (27)) 

angle that system rotates during a quarter cycle at the natural frequency 
(see fig. 8(b)) 

ratios of gondola mass to centrifuge mass (see eqs. (1)) 
ratio of centrifuge mass to space-station mass (see eqs. (1)) 
real part of Laplace variable s 
real part of nondimensional Laplace variable, 
nondimensional time referred to upper critical speed, u> e t 
reference nondimensional time, u> e t 0 (see fig. 8) 
complex modal matrix (see eq. (27)) 

assumed lateral mode shapes for centrifuge 
centrifuge rotational speed relative to space station 
nondimensional rotational speed, 

w e 

rotational speeds corresponding to lower and upper critical speeds, 
respectively (see fig. 7) 

imaginary part of Laplace variable s 
imaginary part of nondimensional Laplace variable, 

w e 

uncoupled gondola-arm natural frequencies Jk^Tin^, respectively 



cog uncoupled rigid centrifuge-support natural frequency, \JK /M c 

^1* > w 2* nondimensional gondola-arm natural frequencies respectively 

co g* nondimensional centrifuge-support natural frequency cog^cog 

co e upper critical speed (frequency) , ^|k(m c + M S ^M S M C 

co . coupled natural frequency of system, i = 1, 2, 3, 4 

n ? i 

co* . nondimensional natural frequency of composite system, i = 1, 2, 3, 4 

n,i 

co 9 uncoupled natural frequency of antisymmetric degree of freedom (see 

ec. (29)) 

Bars over symbols denote Laplace transformed variable. 

Superscript T denotes transposed matrix. 

ANALYSIS 

In this section, a basic mathematical model of a centrifuge-equipped space station 
is described and the governing equations of motion are presented and discussed. The 
equations are nondimensionalized and reduced to an equivalent set of first-order equa- 
tions which are summarized in matrix form. From this matrix form the eigenvalues and 
eigenvectors of the configuration can be easily obtained. 

The present analysis is devoted mainly to an undamped symmetric configuration. 

By using a coordinate transformation and nulling a system mass parameter, the equations 
of motion may be partially uncoupled. From the resulting set of equations, limiting curves 
of the stability boundaries are obtained and interpreted. In addition, the response of the 
configuration is expressed for a given initial state in terms of the eigenvectors and the 
effect of damping and mass asymmetries has been included in the analysis. 

A detailed derivation of the equations of motion is given in appendix A. 

Mathematical Model 

A schematic illustration of the space-station — centrifuge model is given in figure 1. 
The mass of the space station without the centrifuge is denoted by M s . The centrifuge is 
represented by a hub (pivot) mass M^ and two cantilever beam-like arms upon which 
the gondola masses mj and m 2 are attached. The arms of the centrifuge are assumed 
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to be massless with bending stiffnesses of kj and k 2 , respectively. The lengths of the 
respective centrifuge arms are denoted as Lj and L£. 

Assumptions used in the analysis are as follows: 

(1) Motion is planar. 

(2) The space station is kinematically constrained from rotating. 

(3) The centrifuge rotates at a constant angular velocity S2. (This constraint per- 
mits linearization of the equations of motion; however, the angular momentum and energy 
of the system are no longer conserved.) 

(4) Effects of gravitation have been neglected. 

(5) Small motions are assumed. 

(6) The hub mass Mh is isotropically spring mounted and isotropically damped 
relative to the space station with springs of stiffness K and damping C. Damping of 
the centrifuge arms relative to the hub are denoted by Cj and C£. 

Of particular importance are the axial displacements at the ends of the arms which 
result from lateral deflections. (See appendix B.) This effect provides stiffening pro- 
portional to the square of the rotational speed £2 with significant effects on the stability 
of the system. 


Equations of Motion 

Equations of motion have been derived by the method of Lagrange in appendix A for 
a general configuration. Subsequently, the equations are reduced to a simpler form for 
application to the assumed model wherein they become a set of linear differential equa- 
tions with constant coefficients. (See eqs. (A58).) 

Coordinate system .- Motion is referred to a coordinate system X,Y which rotates 
at the same rate £ 2 as the centrifuge. (See fig. 2.) In this rotating frame, the center 
of mass of the space station mass M s (without the centrifuge attached) is located rela- 
tive to a fixed (inertial) frame Xf,Yf by the coordinates x g ,y s and the angular dis- 
placement £2t. As indicated in appendix A, the constant rotational speed constraint 
imposes only an externally applied moment to the system; therefore, the origin of the 
X,Y and Xf,Yf coordinate systems may be conveniently interpreted as being located 
at the center of mass of the composite system. Relative to the center of mass of the 
space station, the displacement of the hub mass M^ is denoted by the coordinates x^y^ 
and finally, the lateral deflections of the gondola masses are denoted as q^ and qg. 

The axial displacements which are directed toward the hub are given by a^ and ^ 
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Reference to equations (A58) of appendix A indicates that the last four equations 
(in coordinates x^, y^, q^, and qg are independent of the coordinates x s and y g . 
Thus, although equations (A58a) and (A58b) are coupled to equations (A58c) to (A58f)), the 
stability of the configuration is determined by the characteristic behavior of the latter 
equations. Therefore, the analysis is concentrated on these relations; that is, equa- 
tions (A58c) to (A58f). 

The Laplace transform is applied to the equations, and to provide a more general 
development, the following nondimensional parameters are introduced into the equations: 


P = 


~\ 




w e 

M c = mj + m 2 + M h 
m l 

m i = m^ 

mo 
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In these relations it should be noted that spring constants kj and k 2 have been sub- 
stituted for the effective lateral stiffness of the beams that are described in equa- 
tions (A17) and the upper critical speed w e has been used as a normalizing factor. 

In order to facilitate an eigenvalue -eigenvector solution, the equations are written in 
first-order form by using the following definitions: 


X 1 

= P x h 

x 2 

= py h 

x 3 

= pq x 

x 4 

= pq 2 

X 5 

= x h 

x 6 

= ^h 

x 7 

= *1 

x 8 

= ^2 , 


( 2 ) 


where the bar symbol (“) has been used to denote a Laplace transformed coordinate and 
p is the nondimensional Laplace operator. Thus using the relationships 


X 1 = P x 5 


x 2 = P x 6 


Xg = PX ? 

x 4 = px g 




( 3 ) 


together with the previous definitions, equations (A58c) to (A58f) may be written in the 
form 

(E + pI)X = 0 (4) 

where E is given for the general case as 
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and I is the identity matrix and 


x = [% } x 2 ,x 3 ,x 4 ,x 5 ,x 6 ,x ? ,x 8 J (6) 

It should be noted that matrix equation (4) is in the form of a nonsymmetric eigenvalue 
problem, the nondimensional Laplace variable p being the eigenvalue parameter. Equa- 
tion (4) represents eight first-order linear differential equations with constant coefficients. 


Undamped symm etric configuration.- The undamped symmetric centrifuge is of 
special interest because by means of a coordinate transformation a partial uncoupling of 
the equations result thereby and provide a deeper insight into the stability characteris- 
tics of the configuration. In addition, as will be shown subsequently, mode shapes of 
undamped configurations may be conveniently displayed. The equations of motion for the 
undamped symmetric configuration may be obtained by specializing the nondimensional 
parameters of equations (1) and the coefficients of the E matrix in accordance with the 
following: 


m 1 = m 2 
L 1 = L 2 


k 


1 





( 7 ) 


h 



o 


One of the four degrees of freedom may now be uncoupled from the analysis by 
combining coordinates. This uncoupling may be done by recognizing that 


p(x 3 + S 4 ) + + ^(LjAi - 1)](, 7 + x 8 ) = 0 

and 

P (X? + Xg) - (xg + x 4 ^ = ° 


( 8 ) 

( 9 ) 


Substituting equation (9) into equation (8) gives the uncoupled second-order equation of 
motion in the new composite coordinate x^ + Xg 



+ A j — 



( 10 ) 


Equation (10) describes an antisymmetric motion of the centrifuge arms about the hub 
which is independent of the other degrees of freedom of the system. Since Lj Aj ^ 1 
for the configuration (as indicated in appendix B) , the quantity in the brackets is always 
positive. Thus, the motion of this composite degree of freedom is always stable for all 
rotational speeds S2* and may be disregarded when stability is considered. 




To complete the uncoupling process, the composite coordinate - Xg is formed 
which remains coupled to the x g and x g degrees of freedom. 

The uncoupled form of equation (4) may now be formulated by using the following 
definitions: 


Z 1 

z 2 

z 3 

h 

z 5 

z 6 

i 7 

z 8 



= pz 5 

= pz 6 
- x 4 = pz 7 


+ x 4 = pig 





+ x 8 


J 


( 11 ) 


and 

Z = |^z 1 ,z 2 ,z 3 ,z 4 ,z 5 ,z 6 ,z 7 ,zg 
Thus, the eigenvalue problem becomes 


( 12 ) 


(E+pl)z = 0 (13) 

where E is the uncoupled form of the previously defined E -matrix for an undamped 
symmetric configuration: 
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0 
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(14) 


Attention will be focused on the three coupled degrees of freedom. 


Characteristic Stability Equation 


The characteristic stability equation for the complete configuration is obtained by 
setting the determinant of the coefficients of either equation (4) or equation (13) to zero, 
that is, 

| E + pi | = 0 (15) 


E + pi = 0 


(16) 


Thus, if a root of equation (15) or (16) has a positive real part, the system is unstable. 

It should be noted that equations (15) and (16) have the same characteristic polyno- 
mial and roots; however, for the undamped symmetric configuration the polynomial of 
equation (10) which is stable may be readily factored from equation (16) and disregarded. 
Thus only the three coupled degrees of freedom need to be considered, that is, the 
unshaded portion of E. By taking the determinant of the remaining 6X6 terms of the 
E -matrix a sixth-order polynomial in terms of z z 2 , z 3 , z 5 , Zg, and z^ results 
of the form 

p 6 + Aj^p 4 + A 2 (j2*)p 2 + A 3 (fi*) = 0 (17) 
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The coefficients Aj, A 2 , and A 3 of the polynomial are functions of the physical param- 
eters describing the system as well as the nondimensional rotational speed £2* . Note 
that the polynomial only contains coefficients of even powers of p as is usual for non- 
rotating systems and roots occur symmetrically about both the real and imaginary axes. 

A three-dimensional plot of the locus of the characteristic roots of a symmetric 
configuration in terms of the nondimensional Laplace variable p and rotational speed 
£2* is shown in figure 3. A modal response of the form Y = Y 0 e^ r has been implied 
where Y is the (complex) response, Y 0 is a complex coefficient, and r is a non- 
dimensional time, r = co e t. The additional dimension £2* has been included to illustrate 
the effect of rotational speed on the roots. Since complex roots occur in conjugate pairs, 
only the positive imaginary roots are shown. In general, roots for the undamped configu- 
ration are stable and have purely oscillatory characteristics; hence, the roots are on the 
jo>*,£ 2 *-plane where the imaginary part of p is o>*, a nondimensionalized natural fre- 
quency of the system. In the analysis for the undamped configuration, it will be advan- 
tageous to examine the locus of the roots on the jo>* £2*-plane, that is, for p = ju>*. It 
should be noted that the natural frequencies co n are frequencies that would be exper- 
ienced in the rotating system as opposed to frequencies observed for the system in a 
fixed reference. 

A typical root as shown at (a) is identified as the lowest of the three (coupled) con- 
ventional natural frequencies which the system has without rotation. If the behavior of 
this frequency is noted while the rotational speed is increased, it goes to zero after which 
it increases and may coalesce with another root before increasing indefinitely. Except 
where the two roots coalesce and the lowest root reduces to zero the roots generally 
increase along the imaginary axis with increasing rotational speed. This trend is due to 
a stiffening effect derived from the centrifugal forces of the rotating system. 

Two unstable regions have been illustrated in figure 3 and the positive portion of 
the surface enclosed by roots which are not purely imaginary is shown shaded. The insta- 
bility which occurs at the lower rotational speed £ 2 * for roots in the region of (b), (c), 
and (d) has no frequency or oscillatory characteristics but a purely divergent behavior 
with real roots as illustrated on the o*, £2* -plane. This behavior occurs between the two 
critical speeds derived from the bidirectional properties of the system and hence has 
been termed a critical speed instability. A second instability identified by roots having 
purely complex values and characterized by an oscillatory divergent behavior may exist 
at still higher rotational speeds. It is analogous to the mechanical instability exhibited 
by helicopter pylon-rotor configurations. (See refs. 2 and 3.) 


* Positive in the sense of having roots with positive real parts. 
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When the undamped nonsymmetrical configuration is analyzed, the characteristic 
polynomial remains even but the equations of motion do not uncouple in the previously 
described manner. Thus, the locus that was ignored in the symmetric analysis must be 
included on the root locus display of figure 3. In fact, this locus may add another oscil- 
latory instability. 

For a configuration where damping has been introduced the characteristic polyno- 
mial is no longer even but, in general, may have coefficients of odd powers in p also. 
Consequently, the root locus will not be symmetric about the imaginary axis. 

The roots of the undamped system have been investigated by using two methods. 

The first method examines the root locus on the -plane of figure 3 and is described 

in the next section. The second approach utilizes a stability criterion which is adapted 
for the even characteristic polynomial to describe unstable regions. The mechanics of 
the adaptation may be seen by noting that a complex or p-plane plot of figure 3 may be 
illustrated by figure 4(a). The same locus in the p2-plane will appear as shown in fig- 
ure 4(b). Therefore, since all stable roots on the p-plane for all anticipated ranges of 
must be on the imaginary axis, all stable roots on the p2 -plane must lie on the nega- 
tive real axis. The Meerov- Fuller criterion (ref. 9) determines when a root of a polyno- 
mial is not negative and real. 

With the addition of damping to the system, the locus plot is no longer confined to 
the jcu*,f2* -plane and neither of the preceding techniques is applicable. However, the con- 
ventional Routh criterion (for example, see ref. 10) is applicable to the configuration and 
has been used to establish unstable regions. 

Method of Uncoupling 

An interesting insight into the nature of the oscillatory instability and the conditions 
under which it occurs can be found by examining an uncoupled form of the stability deter- 
minant of a symmetric undamped configuration on the -plane as has been done in 

reference 2. The uncoupled root loci form limiting curves which approach the coupled 
loci when the rotational speed £1* is large, and these uncoupled loci are shown by the 
heavy solid lines and curves in figure 5. In the figure, centrifuge loci corresponding to 
two different configurations denoted by type I and type H curves have been included and 
the heavy straight lines with slopes of ±1 depict the upper half of the uncoupled hub loci. 

To develop the equations of the loci, it is convenient to refer to the equations of 
motion in second-order form. The matrix formulation for the symmetric configuration 
in terms of Zg, Zg, and z^ is 
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An uncoupled analysis may only consider diagonal terms of a dynamical matrix; 
however, because of the strong coupling of x^ and y^ through the Coriolis accelera- 
tion in this problem, these degrees of freedom must be considered coupled, that is, the 
off-diagonal terms -2ft* p and 2ft* p of the matrix equation (18) must be included. 
Furthermore, the condition that pj — 0 must be impressed on the upper left 2x2 sub- 
matrix. Thus, the uncoupled hub behavior is represented by 

[p 2 + (l - ft* 2 )] 2 + 4ft* 2 p 2 = 0 (19) 

which reduces to 

±(ft*-l), ±(ft* + l) (20) 

o> e 

Sf " COp 

on the jar, ft -plane where p = j and u> n is the natural frequency. These equations 

w e 

correspond to straight lines at ±45° to the ft* -axis of figure 5. The centrifuge behavior 
is represented by a hyperbolic curve corresponding to the diagonal element of the 
Zrj = q-^ - qg degree of freedom. Thus, the equation for this uncoupled locus is 

p2 + r^( L i a i + ^i-0 + rii =o (21) 

The approximate coupled root behavior of the centrifuge interacting with the hub is 
shown dashed at points A, B, and C. When two uncoupled loci cross, the coupled root loci 
may coalesce in a manner similar to that shown by the dashed lines in the neighborhood of 
the intersection of the type I uncoupled centrifuge locus and the ft* - 1 asymptote as 
indicated at point A of figure 5. On the other hand, the root locus may couple without 
coalescing as shown at B, or in the neighborhood of the intersection of the type II uncoupled 
centrifuge locus and the ft* + 1 asymptote as, for example, at point C. An oscillatory 
divergent instability occurs when roots coalesce and no longer remain purely imaginary, 


16 



that is, become complex. No instability results however when the roots couple without 
coalescing. 

It can be seen that when the limiting slope of a hyperbolic centrifuge locus for large 
values of £2* is unity or greater, that is, the uncoupled centrifuge locus resembles that 
shown for type n, no unstable crossings of the hyperbolic curve and the £2* - 1 asymp- 
tote occur. The condition which eliminates the intersection may be found by examining 
equation (21) for large £2, that is, 

p2 + rr^;( L i A i + 2 H ' 1 ) a0 ( 22 ) 

and making the slope equal to or greater than unity. Therefore, the following relation 
must be true: 


L.J A^ + 2p i - 1 

1 - 2m 

Thus, no oscillatory instabilities should exist for 


g 1 




, 2 - L l A l 


(23) 


(24) 


A result similar to expression (24) was established in reference 2. 

In terms of a which is introduced in appendix C (see eq. (C6)) as the proportion 
of the length that the distance of the effective hinge point of the centrifuge arm is from 
the hub, equation (24) becomes 



1 - 2a , 


1 - a. 


(25) 


Thus, for a symmetric configuration having arms behaving as uniform beams, a = 1/6 
and oscillatory instabilities should be prevented when 


§ 0.2 


(26) 


Because of the intrinsic geometrical nature of the parameter a, it is used throughout 
the remainder of the paper in lieu of A. 


Modal Representation 

Since the equations of motion are linear with constant coefficients, it is possible to 
represent the motion of the system by a superposition of the eigenvectors. Of special 
interest is the undamped configuration which exhibits sustained oscillatory motion (imag- 
inary eigenvalues) when it is rotated at a constant angular speed. 
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If the eigenvalues or roots previously found are ascertained to be distinct, eigen- 
vectors which satisfy the equations of motion may readily be found. Since the E -matrix 
as given by equation (14) is real and nonsymmetric, the eigenvectors and eigenvalues of 
equation (13) may be real or may occur in complex conjugate pairs. As in reference 11, 
a general expression for the time history may be given, in terms of nondimensional 
time r, as 

Z(t) = 0e AT 0 -1 Z(0) (27) 

where <fi is the modal matrix and Z(0) is the initial state vector in first-order nota- 
tion, A constitutes a diagonal matrix consisting of the eigenvalues p^ of the system, 
and e Ar is a diagonal matrix having elements e^ T . In general, the modal displace- 
ment will be complex so that conventional time-invariant mode shapes do not exist but 
since eigenvectors occur in complex conjugate pairs and, when superimposed, may yield 
a real response, the time history may be illustrated by real modes which change with time. 
(See ref. 12.) Thus, for the undamped rotating system, its general motion may be dis- 
played with real motions composed of pairs of mode shapes phased in time a quarter of a 
cycle apart. 

When damping is introduced, the motion, in general, cannot be represented with 
real modes. Therefore, no convenient display is available. 

DISCUSSION OF RESULTS 


In this section, results are presented for several space-station — centrifuge config- 
urations in order to illustrate the types of instabilities discussed previously. Physical 
parameters of the system are varied with respect to a reference model. 


Basic Model 

Properties .- The basic model assumed for this study was an undamped symmetric 
configuration defined by the following nondimensional parameters: 


dl = M 2 = °- 1 '' 

jUg = 0.01 

a l = a 2~ 6 
* * 

"i - "a - °-U 


(28) 
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The a values have been assumed to be consistent with those of a uniform cantilevered 

k 1 

beam under the influence of a tip load (eq. (C6)) and a stiffness ratio -=± = = 0.00101 

* K. K. 

has been implied by the co* value. 

Instability reg ions.- A plot of the natural frequencies u>* of the basic model for 
rotational speeds S2* is shown in figure 6. The uncoupled antisymmetric mode (dashed 
curve) which has been neglected in the stability analysis is included for completeness. 

The faint broken lines denote the uncoupled loci of the degrees of freedom considered in 
the stability analysis which were introduced in figure 5. Thus, the basic model exhibits 
a type I centrifuge locus as classified in figure 5, and the frequency plot associated with 
the coupled degrees of freedom resembles that illustrated on the ju>*,fi* -plane of figure 3. 

Inasmuch as the loci corresponding to each degree of freedom of the system must 
be on the -plane for stable operation and the system has four degrees of freedom, 

figure 6 shows two ranges of 12* for which the basic model is unstable. Points of iden- 
tification a, b, d, e, and g of figure 3 are also shown in the figure. The region that is 
encountered at the lowest rotational speed occurs for a zero natural frequency in the 
rotating reference frame and extends from 12* = 0.725 to 12* = 1. (It should be noted 
that the upper boundary is unity because the axes have been normalized to the undamped 
upper critical speed.) The behavior of the degree of freedom in the region is character- 
ized by a nonoscillatory or zero-frequency divergence. The two general displacement 
patterns for the critical speeds which bound the unstable region and occur as a result 
from the bidirectional properties of the configuration have been illustrated in figure 7. 

The various masses are indicated by circles of varying radii. 

A second region where the locus is discontinuous occurs for 2.75 2 12* = 4.22 in 
figure 6. Unlike the first region, this unstable region is bounded by stable loci which 
possess frequency components. It follows that the unstable region as illustrated in fig- 
ure 3, has an oscillatory divergent behavior. The motion, as suggested by the intersec- 
tion of the uncoupled loci of figure 5, is derived from a coupling of the hub and centrifuge 
motions and is analogous to the self-excited mechanical instability experienced with 
helicopters (refs. 2 and 3). Although the usual configuration may, in general, preclude 
this instability, to ensure that its occurrence will not be a factor, the conditions under 
which it may occur have been emphasized in this work. 

Mode shapes .- Further insight into the mechanics of the motion may be gained by 
examining the associated eigenvectors that characterized the motion for a particular value 
of S2*. As previously indicated, the associated eigenvectors are not necessarily real 
and, in general, are complex although they may be combined to represent a real motion. 
The motion may be interpreted as a combination of two conventional modes having com- 
mon frequencies but with a quarter of a cycle phase difference. For example, a typical 
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combined eigenvector or "mode shape” may be graphically interpreted as shown in fig- 
ure 8 where the undeformed schematic representation is shown in figure 8(a) for refer- 
ence. The + symbol, it should be noted, indicates the center of mass of the system and 
all motion is referenced to the rotating X,Y frame as shown. One gondola is marked 
for identification and as previously denoted, the various masses are indicated by circles 
of varying radii. 

Figure 8(b) shows a pair of displacement patterns indicative of a typical response 
in a mode which is derived from combining complex conjugate eigenvectors. The left- 
hand figure has been arbitrarily normalized to a unit gondola-to-hub deflection. The 
right-hand pattern occurs a quarter of a cycle later after the system has rotated counter- 
clockwise through an angle 0. Nondimensional time t = co e t is indicated. Motions at 
other times than indicated are represented by superpositions of these modes. It should 
be pointed out that because of Coriolis acceleration, the hub and space-station traverse 
elliptic trajectories in a mode. 


Before introducing the mode shapes for the three coupled degrees of freedom, the 
response of the uncoupled mode of equation (10) for an undamped symmetrical configura- 
tion is given. The natural frequency of the mode may be obtained from equation (10) as 


* 


o> 


z,2 




+ O 


,*2 


ex i \l/2 


(29) 


where p = ju>* 9 has been assumed and the corresponding mode will have an antisym- 
metric deflection shape as shown in figure 9. 


In order to identify the characteristic behavior of the remaining degrees of freedom, 
especially near an instability, mode shapes for the three coupled degrees of freedom have 
been computed and illustrated in figure 10 for selected (counterclockwise) rotational 
speeds. For each rotational speed £2*, the pairs of modes have been arranged with 
natural frequencies increasing from left to right. At £2* = 0, the classical mode shapes 
which consist of two lateral motions and a longitudinal motion are shown. As the rota- 
tional speed is introduced, centrifugal accelerations couple the longitudinal and lateral 
motions (as shown by the mode shapes for £2* = 0.1 and 0.7). The large angle 6 at 
£2* = 0.7 for the low frequency mode is indicative of the natural frequency approaching 
zero. Actually, at £2* = 0.725 (fig. 6), the frequency co* becomes zero and has the 
nonoscillatory shape as illustrated in figure 7(a). When £2* equals 0.8 and 0.9 (fig. 10), 
the low frequency mode has split into two different modes, the eigenvector has become 
real, and the need to combine complex conjugate modes has been eliminated. The mode 
shown at the left is unstable (divergent) whereas its companion mode on the right is stable 
(convergent). Although the response of these modes changes in an exponential manner, it 
has no oscillatory nature. For £2* = 1.0, the two modes combined to yield the critical 
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speed mode which was indicated in figure 7(b) and when £2* = 1.1, the coupled-mode pair 
relationship which existed at rotational speeds £2* below the unstable region and for the 
other modes in the speed range discussed is reestablished. 

For rotational speeds from £2* = 1.0 to £2* = 2.7, the system is stable and exhib- 
its three modal pairs, but at £2* = 2.7 as indicated by figure 6, the lower frequency 
modes approach a common frequency (and angle 6) and mode shape. For a slight 
increase in rotational speed, these modes couple and one mode which has been shown on 
the left exhibits the unstable behavior found in the oscillatory divergent instability region. 
As indicated, this motion exists for 2.75 ^ £2* ^ 4.22. The mode illustrated in the center, 
for the rotational speeds in this unstable region, is a companion of the unstable mode but 
it is stable. By noting the relative positions of the centrifuge and space station at 
£2* = 2.7 and £2* = 4.3, the figures show that the phasing of the interacting mode pairs 
is reversed as the rotational speed passes through the oscillatory instability region. At 
higher rotational speeds the motion is typical of stable motion. 

Motions which are approached as the rotational speed becomes large (labeled 
£2* — oo) have been included for completeness. The two higher natural frequencies 
approach the rotational speed frequency and the modal motion degenerates to pure rota- 
tion; this condition would be typical of the system if the centrifuge arms were rigid. The 
lower frequency never achieves the rotational speed frequency. Its motion remains oscil- 
latory although the space station motion subsides and only lateral motion persists. 

Effect on rotation .- The effect of rotation on the vibratory motion of the modes is to 
transform lineal oscillations normally experienced in an inertial system into motions 
having both an oscillatory and rotatory nature in the rotating reference frame. In this 
frame the basic motion describes an elliptic path and may traverse the path in a clock- 
wise or counterclockwise sense. This relative motion has been called whirl in the litera- 
ture (see, for example, ref, 2) and may be seen in the mode shape pairs in a reference 
frame that rotates with the centrifuge arms. For instance, when £2* = 0.1, the hub or 
space station in the high-frequency mode traverses a path which is elliptical and opposed 
in direction to the counterclockwise direction of rotation of the centrifuge. 

On the other hand, the intermediate frequency mode exhibits a counterclockwise 
whirl. On examining the whirl of the low-frequency mode, it is shown to vanish as the 
lower critical speed is approached and does not exist in the critical-speed instability 
region. Subsequently, the whirl sense reverses. 

Symmetric Model Study 

Several additional undamped symmetric configurations which represent variations 
from the basic model were investigated in order to ascertain the effects of parameter 
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changes on the stability and vibrational characteristics. Unless otherwise stated, 
parameters of the basic model have been assumed. (See eqs. (28).) 

An operational space station may be designed with heavier gondolas or lighter hubs 
than those prescribed on the reference model. Thus, if the ratio of gondola mass to hub 
mass is made larger than 0.2, the oscillatory instability may be avoided. An analysis of 
a configuration which is the same as the basic model except that p^ is 0.3 has been 
conducted to investigate this case and to compare the behavior of the basic model with a 
configuration without the oscillatory instability. Figure 11 shows the frequency loci for 
the configuration. It should be noted that this configuration has a type II centrifuge locus 
as classified in figure 5. The critical speed range is slightly wider than that for the 
basic model; however, since the uncoupled loci of the centrifuge and hub do not cross as 
shown in region A of figure 5, this configuration is always stable for 12* > 1. The only 
unstable speed range is 0.55 < £2* < 1.0 which is the critical-speed instability region. 

Mode shapes for the configuration are shown in figure 12. They have been included 
to illustrate the interaction of the modes when there is no oscillatory instability. Above 
ft* = 1.0, figure 12 does not display the concurrence of the modal pairs and frequencies 
which occurred for ft* = 2.7 in figure 10 when the oscillatory instability was present. 
Instead, the interaction which occurs for the upper two frequency modes is stable and the 
two lower frequency modes approach hub modes similar to the upper two of figure 10 as 
ft* approaches °°. 

Effect of mass ratios .- A broader picture of the effects of the mass ratio p^ and 
stiffness ratio kj/K on the instabilities of the basic model is shown in figure 13, where 
the Meerov- Fuller criterion adapted for undamped systems has been used. Note that for 
this approach, the centrifuge arms are assumed to have spring stiffnesses = m^co^ 
or k^K = P l w i/ (1 ” P3). Figure 13(a) shows the effect of p^ and kj/K on the 
critical-speed instability region. The ordinate displays the rotational speed normalized 
to the upper critical speed which at ft* = 1 by virtue of the normalization, coincides 
with the upper bound of the unstable region for all p^. Since the lower critical speed 
is the lower bound of the instability region, the system is stable between the critical 
speeds and these regions are indicated for various values of pj. The figure shows that 
although small mass ratios pj may reduce the instability region, small stiffness ratios 
tend to make the region larger. Thus, a larger stable region below the instability regions 
may be achieved for smaller pj and larger k^K where the upper and lower critical 
speeds approach one another. It should be noted that these steps reduce the bidirectional 
or orthotropic characteristic of the centrifuge. 

Figure 13(b) displays the effect of pj and kj/K on the higher speed instability 
region for three values of pj. The system is unstable within the indicated bands. As 
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previously indicated, no oscillatory instability exists for equal to or greater than 

0.2. Figure 13(b) shows that although increasing p^ (below pi = 0.2) tends to delay 
the instability to slightly higher £2* at small kj/K, the regions occur at lower £2* 
for large kj^K and they are wider for all kj/K. Thus larger stable regions above the 
critical-speed instability region may be attained for pj g 0.2 or when p^ is small 
and kj^K is large. 

The effect of pg on the instability regions has not been displayed in the analysis 
since the system instabilities are relatively insensitive to this parameter in the ranges 
of interest. For 0 g pg ^ 0.2, the instability regions changed insignificantly from that 
of the basic model although for pg approaching 1, the instability regions were reduced. 

Effect of hin ge-distance ratio.- Another parameter which affects the stability in a 
manner similar to mass ratio pj is the hinge distance ratio a ^ as given by equa- 
tion (C6). The hinge distance may be changed by altering the stiffness distribution in the 
arms of the centrifuge. The effect of changing a j is illustrated in a plot showing the 
variation of normalized natural frequency co* with rotational speed £2*. Such a dis- 
play is shown in figure 14 for the basic model having an effective hinge point closer to 
the hub, that is, for a\ - 1/12. Figure 14 resembles the plot for the basic model except 
that the oscillatory instability occurs at a lower speed and the critical-speed instability 
region is broader. Figure 15 shows that when = 1/3, that is, the model hinge-point 
distance has been doubled, the natural-frequency variations show no oscillatory instability, 
and the critical-speed instability region is less. Indeed, as shown in the analysis, by a 
proper choice of pj or otj, the oscillatory instability may be eliminated. The outboard 
choice of hinge point, when it is feasible, appears to offer improved stable operating 
ranges. 


Unsymmetrical Model Studies 

By altering the relationships between pj and pg, aj and otg, or co!£ and coi^, 
models which lack mass symmetry about the hub may be obtained. It should be noted, 
however, that mass balance, as reflected by equation (A57), was retained. In this study 
only mass asymmetries are considered which vary the natural frequencies as 



co 


* 
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This relationship may be derived by assuming that both centrifuge arms are uniform 
beams with equal flexural rigidities. 

Two unsymmetrical configurations were studied. The stability characteristics of 
the first configuration where the effective hinge distance of one arm a j has been halved 
from that of the other arm are shown in figure 16. Because of this change, a second 


/ 
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instability region at S2* = 1.5 having an oscillatory divergent nature is introduced. 
Furthermore, figure 16 shows that the uncoupled mode no longer exists but becomes 
coupled into the system. 

The second configuration may resemble a system having a centrifuge with a single 
gondola. For this model, pj was halved (and doubled) and the natural frequencies 
are given in figure 17. The dynamic characteristics appear to be similar to those of 
figure 16 although in figure 17 the oscillatory instability regions are not as wide. Indeed 
the lower speed region is extremely narrow. 

Thus, a general configuration may be conceived of having not only one but possibly 
two oscillatory-divergent instabilities whereby the condition of symmetry removes one 
of them. Since most actual configurations may have slight asymmetries, more than one 
oscillatory divergent instability region may be anticipated. 

These studies on an undamped system indicate that unstable regions may occur for 
some rotating mechanisms such as the space-station centrifuge. Although the critical- 
speed instability might be anticipated, there is a second type of instability which could 
represent a dangerous operational condition and must be considered. Even though the 
instability always occurred at rotational speeds above the critical speeds, it has been 
shown that such an instability need not exist for an undamped system. In the next section 
an attempt is made to ascertain some effects of damping on the unstable regions. 


Introduction of Damping 

In an effort to ascertain trends that may be expected when an energy dissipation 
mechanism is added to the symmetric mathematical model, viscous dashpots were added 
between the gondola masses and the arms and between the hub and support structure 
(see eq. (A19)). The motivation for considering damping effects stems from the fact that 
damping can exert a destabilizing influence on systems which are nonconservative with- 
out damping (ref. 13). Thus, unstable regions for a series of damped configurations of 
the basic model were determined by using the general equations (5) and the Routh sta- 
bility criterion. Hub damping values of £ = 10“®, 0.01, 0.1, and 0.2 and gondola 
damping values of = f? 2 = 0.01, 0.1, and 0.5 were inserted in the model, and stability 
was determined for a series of and kq/K ratios. Values of the critical damping 
ratio c \/C may be obtained from the figure by using the following equation: 


fi-iivES 

C ?! K 


(30) 


The studies conducted with equal nonzero gondola damping values = ? 2 ) showed 
that gondola damping had no effect on the critical -speed instability. The hub damping 
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ratio £, however, had considerable influence on this instability. For example, with 
£ = = 0.01, the frequency loci in the critical-speed region differed little from that 

shown in figure 13(a). However, increasing the hub damping ratio £ resulted in con- 
siderable reduction in the purely divergent instability regions as shown in figure 18. As 
may be seen in this figure, this unstable region may be eliminated for proper combina- 
tions of gondola mass and stiffness ratio when adequate hub damping is present. 

Furthermore, the results of the analysis indicate that both types of damping influ- 
ence the oscillatory divergent instability region. For large, as when £ = 10“® 

and = 0.01 or 0.5, the stable regions above the critical-speed instability region were 
very small for all configurations with 

0.01 g jUj g 0.49 

Equation (26) established a critical value of p^ above which the oscillatory instability 
ceased to exist for undamped configurations. However, since p^ for a symmetric con- 
figuration must be less than 0.5, this inequality indicates that for a damped configuration 
where is large, no such limit exists. 

Increasing results in an increased or decreased oscillatory divergent insta- 
bility region. For example, consider the basic model with £ = 0.01. When = 0.01, 
the unstable regions are reduced only slightly from the model with ( = (j a ^ = 0, but 
when = ?2 = 0*1 > closed oscillatory unstable regions were obtained as shown in fig- 
ure 19. The unstable region when = 0,01 has become insignificant in the ranges of 
stiffnesses and rotational speeds considered, and more configurations are shown to be 
stable. Although not shown, further reductions in unstable regions were found when the 
hub damping ratio £ was increased to 0.1 and 0.2. Furthermore, the unstable bands of 
the undamped model shown in figure 13(b) for p^ = 0.05 and 0.1 do not extend to higher 
stiffness ratios in figure 19, but generally encompass a larger range of rotational speeds 
for lower stiffness ratios. Figure 19 further shows that when p^ = 0.1, there is a slight 
reduction in the £2* range of the unstable region as kq/K approaches zero. This 
trend becomes much more predominant (fig. 20) when is increased to 0.5; as a 
result, there is an extension of the unstable regions over a wide range of rotational 
speeds and a narrow range of stiffness ratios. Of special interest is the appearance of 
the unstable region for = 0.2 which did not exist for the model with less damping. 
Furthermore, this region represents a region which is isolated to a finite range of stiff- 
ness ratios. An additional value of pj_, that is, p-^ = 0.15, has been included in the 
figure to indicate the increase in size of the unstable region which can occur as the 
(pi = 0.2) critical value of p ^ for the undamped model is approached. 



RESUME 


A linear analysis has been conducted to determine the vibrational characteristics 
and stability boundaries of a centrifuge — space-station configuration. The mathematical 
model used in the analysis describes a centrifuge represented by two gondola masses and 
a hub mass, rotating at constant speed relative to a free but nonrotating space station. 
Massless beams which were cantilevered from the hub constituted the centrifuge arms 
and were assumed to support the centrifuge gondolas. The hub was attached to the space 
station by an isotropic spring. Planar motion was assumed and gravitational effects 
were neglected. Viscous damping was introduced into the centrifuge arms and support 
structure to study damping trends. Of special importance for stability considerations 
has been the inclusion of an axial component of motion for the gondolas resulting from 
lateral motion of the centrifuge arms. This secondary effect of the motion causes the 
natural frequency of the centrifuge arms to increase with rotational speed and is the 
source of a so-called centrifugal stiffening effect. 

Equations of motion were derived by the method of Lagrange subject to two kinemat- 
ical constraints: (1) the space station was constrained from rotating and (2) the centri- 

fuge was constrained to rotate uniformly. The equations were referred to a coordinate 
frame which rotated uniformly with the centrifuge and which had its origin located at the 
center of mass of the composite system. 

The resulting second-order differential equations were linearized by neglecting 
small nonlinear terms. These equations were nondimensionalized and reduced to an 
equivalent set of first-order equations and summarized in matrix form. Two analytical 
techniques were applied to determine the stability of the system. Stability of undamped 
configurations was studied by examining the natural frequencies at various rotational 
speeds. Instabilities using this technique were indicated by an absence of frequencies 
from the imaginary axis of a complex-plane plot of the roots. On the other hand, the 
effect of parameter changes on stability was demonstrated better by using stability cri- 
teria consisting of either an adaptation of the Meerov- Fuller criterion (for undamped 
systems) or the Routh criterion (for damped systems) . 

As in the case of the restrained helicopter, the study has shown that a centrifuge — 
space-station configuration can exhibit two types of instabilities. The instability which 
occurs at lower rotational speeds has been called a critical speed instability and has a 
nonoscillatory divergent characteristic in a reference frame rotating with the centrifuge. 
The higher speed instability was oscillatory divergent. 

For symmetric configurations without damping, the critical speed instability, which 
was attributed to the bidirectional properties of the centrifuge, was found to exist for all 
configurations. By adjusting parameters to make the centrifuge characteristics less 
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directional, either by concentrating the centrifuge mass nearer the hub or by making the 
centrifuge arms stiff compared with the support stiffness, the region of instability could 
be reduced. Limited analyses indicated that this instability could also be alleviated by 
increasing the arm stiffness near the hub. 

The oscillatory divergent instability may be attributed to the coupled motions of the 
centrifuge and hub degrees of freedom. Only a centrifuge having a ratio of gondola mass 
to hub mass less than a critical value exhibited this behavior and the range of unstable 
rotational speeds was largest for configurations near the critical parameter value. 

Larger critical values resulted when the region of the centrifuge arm nearest the hub 
was flexible. A second range of rotational speeds which exhibited the oscillatory insta- 
bility existed for mass-balanced asymmetric configurations. 

When viscous damping was introduced into the centrifuge arms and support struc- 
ture of the model, both regions of instability were affected. Damping in the support 
structure reduced the critical speed instability and damping in the arms tended to stabil- 
ize certain systems operating at rotational speeds near an oscillatory instability. How- 
ever, excessive damping in the arms had a destabilizing effect at higher rotational speeds. 

Only configurations with a centrifuge constituting only a small fraction of the total 
space station mass were considered. Within this restriction, the dynamic behavior was 
found to be insensitive to variations of these masses. 

Since the matrix form of the undamped model was nonsymmetric, the solutions 
possessed eigenvalues and eigenvectors which were either real or occurred in complex- 
conjugate pairs. A display of such modes of the model were presented. When the roots 
were complex (conjugate), the response appears as a dual mode exhibiting two interacting 
conventional modes. The complete set of eigenvectors was used in a general matrix 
expression for the time history of the resulting real motions. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., March 9, 1970. 


APPENDIX A 


EQUATIONS OF MOTION 

In this appendix, planar equations of motion for a general mathematical model of a 
centrifuge -equipped space station are derived and summarized in a reduced form. The 
equations are derived by the method of Lagrange subject to certain kinematical 
constraints. 

Symbols 

a j(p,t),a 2 (p,t) axial displacements of centrifuge arms due to lateral deflections 
C effective hub support viscous damping coefficient 

D dissipation function 

EIj(p) ,El 2 (p) flexural rigidity of centrifuge arms 

I 0 polar mass moment of inertia of space station (less centrifuge about its 

center of mass) 

Ij^Ig lateral modal mass unbalance about undeformed position (see eq. (A26)) 

i ,j unit vectors along the rotating coordinate axes X and Y, respectively 

K effective hub support stiffness coefficient 

L^,L 2 lengths of centrifuge arms 

Mj,M 2 total mass of centrifuge arms (including hub) (see eq. (A22)) 

M c total mass of centrifuge, + M 2 + 

Mft mass of hub 

M s mass of space station 

m^(p),m 2 (p) distributed mass of centrifuge arms 
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mj,m 2 

m 1 ; m 2 

Ni,N 2 

N xs ,Ny S 


Qx 

q x (t) ,q 2 (t) 


mass of gondolas 

generalized masses (see eq. (A16)) 

nonlinear terms (see eqs. (A36) and (A38)) 

nonlinear terms (see eqs. (A30) and (A33)) 

generalized force associated with x generalized coordinate 

generalized coordinates associated with lateral mode shapes 0^(p) 
and 0 2 (p) 


Rl,R 2 

R s 

S 1’ S 2 

T 


t 


position vectors locating arbitrary points on centrifuge 
position vector of space station center of mass (less centrifuge) 
static mass unbalances of centrifuge about hub (see eqs. (A26)) 
total kinetic energy of configuration 
time 


u 1 (p,t) } u 2 (p,t) lateral deflection of centrifuge arms relative to undeformed position 

V total potential energy of configuration 

X,Y coordinate axes which rotate with centrifuge (see fig. 2) 

Xf,Yf space-fixed (inertial) axes (see fig. 2) 

x arbitrary generalized coordinate (see eq. (A20)) 

Xf( ) ,yj( ) coordinates along fixed (inertial) X^Yj axes 

x h ,y h generalized coordinates of centrifuge hub relative to space station 

(see fig. 2) 

x s ,y s generalized coordinates of space station (less centrifuge) (see fig. 2) 
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j3 generalized coordinate of rotation of centrifuge relative to space station 

y ■= 6 + P 

A^(p),A 2 (p) axial modal displacement factors of centrifuge arms (see eqs. (A6)) 

5 (p) Dirac delta function 

equivalent viscous damping ratios for centrifuge arms 
6 generalized coordinate of rotation of space station relative to inertial space 

X dummy variable 

p running coordinate along undeformed centrifuge arms originating at hub 

0j(p),^2(p) assumed lateral mode shapes for centrifuge (see eqs. (Al)) 

£2 centrifuge rotational speed relative to space station 

o>pU >2 uncoupled gondola-arm natural frequencies 

Dots over symbols indicate time derivatives. 

Mathematical Model 

The space- station centrifuge model shown in figures 1 and 2 consists of a space- 
station mass denoted by M s whose mass center is located relative to an inertial ori- 
gin O by the generalized coordinates x g ,y g measured along rotating axes X,Y with 
respect to the inertial axes x^y^. The hub (center of rotation) of the centrifuge denoted 
by mass M^ is attached to the space station by linear springs and linear viscous dash- 
pots, and the support structure is represented in a linear isotropic manner. The hub is 
located relative to the mass center of the space station by the generalized coordinates 
x h’^h" ^he arms of the centrifuge are represented by two massless beams cantilevered 
at the hub and supporting the gondola masses m^ and mg. The lateral deflections of 
the arms are approximated in terms of distributed coordinates as follows: 

Uj(p, t) = <£ ± (p) q x (tfl 
u 2 (p,t) = 4> 2 (p) q 2 (t)J 
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where q^(t) and q£(t) are the generalized coordinates associated with the mode shapes. 
The quantities <p-y(p) and $ 2 (p), respectively, of equations (Al) may be interpreted as 
approximating the lateral deflections of each centrifuge arm by a single term (the funda- 
mental) in a modal series expansion. 

Motion is allowed only in the plane of the centrifuge and the effects of gravitation 
are neglected. Initially, the formulation allows a rotational degree of freedom 6 about 
the space-station center of mass and a rotational degree of freedom /3 representing an 
angular displacement of the centrifuge about its hub relative to the space station. Later, 
to make the equations amenable to analysis, the motion will be constrained in accordance 
with the following equations: 

e(t) = o\ 

> < A2 ) 

/3(t) = fit) 

where is a constant representing the steady spin rate of the centrifuge. The physical 
interpretation of equations (A2) is discussed in a subsequent section. 


Energy Formulation 

Vector relations .- In terms of the coordinate system (see figs. 1 and 2) which 
rotates with the centrifuge, the position vector R s locating the center of gravity of the 
space-station mass M s may be written as follows: 

R s = Xs 1 + y s J r ( A3 ) 

Similarly, the position vectors and R 2 locating an arbitrary point on the left and 

right arms (beams) of the centrifuge may be expressed as 


(A4) 

In equations (A4) p represents a running coordinate originating at the hub where it is 
positive in the sense of i for the right centrifuge arm and positive in the sense of -i 
for the left arm. The coordinates u^(p,t) and Ug(p ,t) represent the lateral deflections 
of the centrifuge arms and x^y^ locates the hub (center of rotation) of the centrifuge 
relative to the space station. 

It should be noted that small axial displacements a^ and ag have been included 
in the analysis which are representative of the inward displacement experienced by a 
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beam deflected laterally. An expression of this displacement for small lateral deflections 
may be expressed as follows (see p. 252 of ref. 14): 


a l (p ’ t)= lI 0 P 

^(pA) = |J o 


du 1 (X,t) _ 

2 "" 
dX 

dX j 

du 2 (X,t) 

2 

dX 

J 

dX 


(A5) 


where X is a dummy variable replacing p. By using the assumption of equation (Al) , 
this inward displacement (or shortening) may be expressed as 


where 


a 1 (p,t) = | A x (p) q^(t)j 
a 2 (p,t) = | A 2 (p) q 2 (t) 

" J 



(A6) 


(A7) 


It should be noted that although a^(p,t) and a 2 (p,t) may be small, they provide a cen- 
trifugal stiffening term which has a significant effect on stability. 

Equations (A3) and (A4) can be differentiated to give the velocity vectors which are 
required to determine the kinetic energy. In performing this differentiation, it must be 
noted that the unit vectors i ,j are time dependent by virtue of the rotating reference 
frame. If y is defined as 


y = o + p 


(A8) 


then it follows that 


di *“ 

dT = w 








(A9) 


By using equations (A9) and taking the derivative of equation (A3) , the velocity R s of 
the space station is given by 
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R S = (is - ry$ + (y s + y*s)j r 


(A10) 


and from equations (A4) with the use of equations (A9) the velocity of a point on the cen- 
trifuge becomes 

R i = A s + x h - a^t) - r[y s + y h + u !(p,t)jj I 


or 


(All) 


+ |X S + y h + UjO), t) + y[x s + x h + p - a 1 (p,t)jJ j 7 

t 2 = + x h + a 2 (p,t) - r[y s + y h - u 2 (p,t)j| i 

+ (js + y h ' tyPfi + y[ x s + x h - p + a 2 (p,t)jj j" 


R 1 = jjs + x h ' A l (p) q l q l ' ^[ y s + y h + ^l (p) q l]j> i 
+ {*s + y h + q l + ^[ x s + x h + p ‘ | A l(p) q l_ j 1 

r 2 = + x h + A 2 (p) q 2 q 2 - yjy g + y h - <*> 2 (p) q 2 ]Ji 

+ { y s + y h " *2^) q 2 + rj*s + x h - p + I a 2 (p) <lfj} 1_ 


when equations (Al) and (A6) are inserted. By using these relations (eqs. (All) and (A12)) 
the complete kinetic energy of the configuration may readily be formulated. 

Kinetic energy .- The kinetic energy of the configuration consists of the kinetic 
energy of the centrifuge and space station and may be written as 


(A12) 


T = \ M s x 2 + y 2 s + r 2 (xg + yg) + 2y(x s y g - y g x s ) 


+ | l o 0 2 + 1 R 1 • R 1 mi(p) dp + | J 2 r 2 • r 2 m 2(p) d P (A13) 
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where I Q is the polar mass moment of inertia of the space station about its center of 
mass. Thus, when equations (A12) are substituted into equation (A13), the complete 
kinetic energy expression results. 

Potential energy .- The potential energy V for the configuration consists of the 
strain energy due to the displacement of the centrifuge (hub) relative to the space station 
and the strain energy associated with lateral deflection (bending) of the centrifuge arms. 
In deriving an expression for the former, it may be convenient to reference the hub dis- 
placement to a nonrotating frame X,Y fixed in the space station and to use the equations 
for transformating the X,Y axes to the rotating axes, that is, 

x. = x^ cos y - y, sin y 1 

h h } (A14) 

y h = x h sin y -i- y h cos yj 

to obtain the potential energy function for the configuration as given by 

v = | K ( x h + yj) + | “i w i ^(t) + | ffi 2 w 2 q 2 (t) ( A15 ) 


In equation (A15), the constant K is the equivalent isotropic spring constant for the 
centrifuge hub support structure. The other quantities are related to the mode shapes 
used in equations (Al) to describe the lateral deflections of the centrifuge arms. Specif- 
ically, the relations 


m l 

m 2 



(A16) 


are the respective generalized masses associated with each mode shape and 



are the squares of the natural frequencies. (See ref. 4.) 


(A17) 
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Dissipation function .- In a manner similar to that used to derive the potential 
energy, the dissipation function D for the configuration may be initially expressed as 

D = | C (^h + y£) + + “^2^2 (A18) 


Then by using equations (A14) the dissipation function may be written in terms of quanti- 
ties related to the rotating axes (see ref. 3): 


D = | c [xh + Yh + r 2 (x£ + y 2 ) + 2y(x h y h - Vh)] + m^c^q 2 + m 2 C 2 a> 2 q 2 (A19) 


where C is the viscous damping coefficient associated with the hub support structure 
and and £ 2 are the equivalent viscous damping ratios associated with the respec- 
tive arms of the centrifuge. 

It should be noted that because of the homogeneous nature of the hub damping 
assumed for the model, the dissipation is expressed in terms of the total velocity of the 
hub relative to the space station which involves the rotational rate y. 


Kinematical Constraints 

Thus far the space-station — centrifuge configuration has been allowed eight degrees 
of freedom, namely: x s , y s , X^, Y^, q^, q 2 , 0, and /3. It is at this stage of the 
development that the kinematical constraints given by equations (A2) are introduced into 
the formulation. This procedure may be accomplished by either explicitly defining the 
generalized forces (external) which maintain the constraints (Newtonian approach) or by 
implicitly accounting for the effect of these generalized forces by substituting the con- 
straint conditions (eqs. (A2)) directly into the equations for the kinetic and potential ener- 
gies and the dissipation function (method of analytical mechanics, ref. 15). The founda- 
tion for this latter approach which is used in the report utilizes Hamilton's principle and 
the calculus of variations. Thus, only the kinematical constraints need to be specified and 
not the generalized forces that maintain these constraints. It should be noted that with 
this approach the configuration possesses only six degrees of freedom, namely, x s , y s , 

x h> Y h> <*1’ ““d <*2* 

If the generalized forces which maintain the given kinematical constraints are of 
interest, they can be obtained by the method of Lagrange multipliers. (See ref. 15.) For 
the present situation this approach readily shows that the generalized force that maintains 
the constraint j3 = fit is a moment acting at the hub of the centrifuge and similarly that 
the constraint 0 = 0 is maintained by a moment acting on the space station. 
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Finally, since only external moments act on the configuration, linear momentum is 
conserved. The mass center of the configuration therefore remains fixed in space or 
translates with a constant velocity. Hence, the origin O of the space-fixed axes Xf,Yf 
(see fig. 2) may be conveniently interpreted as being located at the center of mass of the 
configuration. 


Lagrange’s Equations of Motion 

The equations of motion for the space-station — centrifuge configuration are derived 
by the method of Lagrange by substituting equations (A13), (A15), (A19), and Q x = 0 into 


±(&T\ .8T av ®_ Q 

dt\8x/ dx.dx.dk S 


(A20) 


where x represents an arbitrary generalized coordinate selected from the set (x s , y s , 
x h, y^, q^, and q 2 ). The technique is illustrated for the choice x * x s . From 
equation (A13), 

• • 

1*2 8Rc 


f L l£ 

c s « “Ag 

and using the results of equations (A12) and defining the following integrals: 


rVT / \ r^l -i- OKl P J-tO C7XVO 

-f - = M s (x s - yy ) + Rj • — ] - rn^p) dp + 'r 2 '^ m 2 (p) dp (A21) 

tes K ' J 0 9x c J o ax 0 


Mj = 

r L i 1 

m i(p) dp 
J 0 1 

m 2 = 1 

p L 2 

m 2 (p) dp 

2 J 


(A22) 


equation (A21) reduces to 
9T 


9x c 


M s (x s - yy s ) + Mijig + x h - y( y s + y h jj - q^j J 1 Aj(p) m 1 (p) dp 

pLi ~ L 2 

- J Q 0j(p) « 1 1 (p) dp + M 2 ji s + x h - y(y s + y h jj + q 2 q 2 j A 2 (p) m 2 (p) dp 

p^2 

+ yq 2 J <M p) m 2^ p) dp 


(A23) 
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Differentiating equation (A23) with respect to time produces 

- w B - yy s ) + Mips + X h - r( y s + y h ) - r(y s + y h )J 

- (q? + q^) J o A l(p) m^p) dp ■ yqi J 0 1 01 (p) m i(p) d P 

- ^1 J 0 1 0i(p) m i(p) d P + m 2 ^s + x h - r(y s + y h ) - y(y s + y h Jj 

/ n \ pijg ^^2 

+ ( q 2 + q 2^2/J Q A 2^P) m 2<P) d P + ^2 J 0 02^ m 2^ dp 

p L 2 

+ yq 2 J 0 0 2 (p) m 2(p) d P 

Next, differentiating the kinetic energy with respect to x s gives 


_8T 

ax c 


= M s (y 2 x s + yy s ) + Mjy^g + y h + y(x s + x h JJ + 1 1 yq 1 + S^ 2 - | y 2 q 2 
x J Q A l(p) 0i(p) m i(P) d P + M 2 0[y s + Vh + 0( x s + x h)] “ *2^2 


- S 2 y 2 + | y 2 q 2 J A 2 (p) 0 2 (p) m 2 (p) dp 


where 


~ r L i r L i "" 

1 1 = Jq 01<P) m l<p) d P S 1 = J 0 P m l(p) d P 

^ pL*2 

1 2 = j 02<P) m 2<P) d P S 2 = J Pm 2 (p) dp 




(A24) 


(A25) 


(A26) 
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From the potential energy function V as given by equation (A15) there is a contribution 
of 



and from the dissipation function D (eq. (A19)), there is a contribution of 


(A27) 


= 0 (A28) 

9X S 

Combining equations (A24), (A25), (A27), and (A28) in accordance with equation (A20) gives 
the Lagrangian equation of motion for the x s degree of freedom (generalized coordinate) 
as follows: 

For the x s equation: 


M. 


; (x s - 20y s - 0 2 x g ) + (Mi + M 2 ) 


x s + x h - 20(y g + 


y h )-n 2 (x s+ x h j 


- (Sj - S 2 )0 2 - 2fil 1 q 1 + 2S2I 2 q 2 = -N xs 


(A29) 


where N xs is composed of nonlinear terms as follows: 

, 0 1a 1 (p) m 2 (P) dp + (q| + ^2)^ 


N xs = "(qJ + 1 A 1 (p) m i(p) d P + (q 2 + q 2 <i 2 ) J 2 A 2^ m 2^ dp 


2 ° 2q l I 0 1 Al ^ dp " I ° 2q2 I 0 2 A2 ^ ^ 2 ^ m2 ^ dP ^ A3 °^ 


In a similar manner, the equations of motion for the remaining degrees of freedom 
may be obtained and are: 

For the x^ equation: 


( M 1 + M 2)p s + x h - 2fi(y s + y h ) - fi 2 (x s + xj 


- - S 2 ^0 2 - 2S2Ijqj 


+ 2l2I 2 q 2 + C(x h - fiy h ) + Kx h = -N xs 
For the y g equation: 


(A31) 


M s (y s + 2 Qh s - 0 2 y g ) + (Mj + M 2 )|^ s + y h + 20 (x s + x h ) - 0 2 (y g + y h ) 

+ 1^ - 0 2 q^ - I 2 (q 2 - 0 2 q 2 ) = -N yg (A32) 
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where 


Nyg — -2fiq^qj Aj(p) m^(p) dp + 2J2q 2 q 2 A 2 (p) *** 2 (p) d P (A33) 

For the equation: 

(M x + M 2 )^y s + y h + 2S2 (x s + x h ) - n 2 (y s + y h ) + - n^i) - i 2 (q 2 " ° 2q 2 ) 


c (y h + « x h ) + K yh = ' N ys 

For the qi equation: 

m^i + 2in 1 ? 1 a) 1 q 1 + 


q l 


2 o o P^l 

jOJj - min' 5 + SI pAj(p) mj(p) dp 

+ ii[y s + ^h + 2 °( k s + k h) - n2 (ys + y^] - q i[ k s + * h - 2n (y s + *h) 

- n 2 (x g + x h) ] 1 Ai(p) m^p) dp = -Ni 


where 


N i = ( q i q i + q i q i ■ I n2q ?) J 0 1 A i^ m i0>) d p 


For the q 2 equation: 


m 2 q 2 + 2m 2 ? 2 0J 2 q 2 +|m 2 a)| - m 2 n 2 + n 2 pA 2 (p) m 2 (p) dp q 2 

- ^[y s + ^h + H k s + k h) - fi2 (y s + y h) + q 2[*s + x h 

- 2n(y g + y h ) + « 2 (x s + x h jj J 2 A 2 (p) m 2 (p) dp = ~N 2 


(A34) 


(A 3 5) 


(A36) 


(A37) 


where 

N 2 = ^q 2 4 2 + q 2 q 2 - | ^ 2q 2 ) J Q 2 a 2 (p) m 2 (p) d P ( A38 ) 

It should be noted that residual nonlinear terms in equations (A3 5) and (A37) are under- 
lined and have purposely not been included in Ni and N 2 . It will be shown later that 
these residual nonlinear terms combine to produce a linear term proportional to Sr, 
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Reduction of the Equations 

In reducing the equations of motion to a form more amenable to a linear stability 
analysis, the nonlinear terms N xs , N ys , Nj_, and N 2 are assumed to be small. Fur- 
thermore, the centrifuge is assumed to consist of two massless arms with bending rigidi- 
ties (EI)j and (El) 2 , a hub mass M^, and two gondola masses and m 2 as shown 

in figures 1 and 2. The mass distributions reflecting these assumptions are given by 

m x (p) = | M h 6(p) + m 1 6(p - L x ) 

m 2 (p) = \ Mh 6(p) + m 2 6(p - L 2 ) 

where the Dirac delta function 5 has been employed in the definition. In addition, the 
mode shapes are assumed to be normalized so that 



0l(Ll) = 0 2 (L2) = 1 

By use of equations (A39) and (A40), it now appears that 
M 1 = \ M h + m l 
M 2 = | M h + m 2 




(A40) 


Sj = m^L^ 

s 2 = m 2 L 2 
h = m i 


I 2 = m 2 


nij = mj 


m 2 = m 2 


nijl 

A x (p) mj(p) dp = mj A^LjJ 

C ^2 

J o A 2 (p) m 2 (p) dp = m 2 A 2 (L 2 ) 


(Equations continued on next page) 
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f* L l 

J 0 P A l(P) m i (p) d P = m i L l A l( L l) 

p L 2 

J o pA 2 (p) m 2 (p) dp = m 2 L 2 A 2 (L 2 ) 
J 0 <Mp) a 1<P) m l^p) d P = m l A l( L l) 

p-^2 

J o 0 2 (p) A 2 (p) m 2 (p) dp = m 2 A 2 (L 2 ) 


By using equations (A41) , the equations of motion may be summarized 
For the x s equation: 

M s (x s - 2J2y g - fl 2 x s ) + (M h + m 1 + m 2 )|x s + x h - 2S2(y g + y ^ - Sl 2 (x s 

- ^mjLj - m 2 L 2 ^fi 2 - 2m^J2qj + 2m 2 S2q 2 = 0 
For the x^ equation: 

( M h + m l + m 2)[ x s + x h ' + K) - fi2 ( x s + x h)] - ( m l L l “ m 2 L 

- 2mjS2q^ + 2m 2 0q 2 + C^x h - fly^ + Kx h = 0 
For the y s equation: 

M s(y s + 2 S2i s - fi 2 y g ) + (M h + m 1 + m 2 ^y g + y h + 2fi(x g + x h ) - fi 2 (y s 

+ - m 2^2 " fi2q 2) = 0 
For the y^ equation: 

( M h + m 1 + m 2 )^ g + y h + 2S2(x s + x h ) - ft 2 (y g + y h J + m^ - S2 

- m 2^ q 2 “ ° 2q 2) + c (yh + ^h) + = 0 


(A41) 


as follows: 

‘ + x h)] 

(A42) 

2 )n 2 

(A43) 


S + 5-h)] 

(A44) 

2 4 

(A45) 
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For the q^ equation: 

nyij + 2m 1 p 1 c t > 1 q 1 + jm^a; 2 - fi 2 ) + mjLj A^L-^J^jqj + m.^ + y h 

+ 2fi (x s + x h ) - fi 2 (y s + y h )] - m 1 A i( L i)|\ + x h 

- 2S2(y g + - fi 2 (x s + XjJ^ = 0 (A46) 

For the q 2 equation: 


m 2 q 2 + 2m 2 p 2 co 2 q 2 + m 2 (o> 2 - ft 2 ) + m 2 L 2 A^L^ft 2 ^ - + y h + 2ft(x s + x h ) 

“ ° 2 (y S + y h)] + m 2 A 2( L 2) f s + % - 2ft(y g + y h ) + ft 2 (x g + x h )Jq 2 = 0 (A47) 


The acceleration coupling found in these equations may be removed by combining 
equations (A42) to (A47); for example, subtracting equation (A43) from (A42) and defining 

M c = Mfo + m^ + m 2 (A48) 

gives the following for the x s equation: 

M s (x g - 2fty g - ft 2 x s ) - c(x h - fty h ) - Kx h = 0 (A49) 

In a similar manner, subtracting equation (A45) from (A44) gives for the y s equation: 

M s (y s + 2ftx g - ft 2 y g ) - c(y h + fix h ) - Ky h = 0 (A50) 

In order to simplify the equation, equations (A43) and (A45) are used to eliminate the 
terms + % - 2ft(y g + y h ) - ft 2 (x s + x h jj and ^ g + y h + 2 ft(x s + x^j - ft 2 (y s + y h jj, 
respectively, from equation (A46) to obtain 

m^Mh + m^q.^ + 2m^M c £ 1 a> 1 q 1 + - m^M^ + m^ft 2 + mjM c Lj 

- m 1 AjfLj) ( m i L l - m 2 L 2 ) J2 ^] <1 l + m l m 2(^2 " ° 2q 2) " m l C ( y h + ^h) “ m l Ky h = 0 

(A51) 

In a similar manner, equation (A47) may be put in the following form: 


42 



APPENDIX A 


i 


m 2(^h m l) ( l2 2m2M c ^2 Ct) 2^2 |^^2^^^2 ~ **■ ^2^c^ J 2 ^2 ^2)^^ 

- m 2 A 2 (h^j(m 1 L 1 - nigL^^qg + mjm^ - 0 2 q^)+ m 2 C(y h + J2x h ) + m 2 Ky h = 0 

(A52) 

It should be noted that in performing the operations to obtain equations (A51) and (A52) , 
the residual nonlinear terms have been retained with linear terms proportional to fi 2 . 

Solving for q 2 in equation (A52) and substituting the resulting expression into 
equation (A51) gives the equation: 

mjMjjMcqj + 2m 1 M c ^M h + m^jWjqj + + m^w 2 - mjMjjMc^ 2 


+ m 1 M c (M h + m^Lj A^L^tt 2 - m^Mh + m^A^L^mjLj - m^fi^jqj 
- 2mjm 2 M c ? 2 ^2^2 " ^nigMc^ + m^m^cL^ A 2 (l 2 )£2 2 

l m 2 A 2( L 2)( m l L l “ m 2 L 2) fi ^2 " m l M c C (y h + ^h) - m l M c K y h = 0 


- m 


(A53) 


Solving for q^ from equation (A51) and substituting the resulting expression into equa- 
tion (A52) gives the q 2 equation: 

m 2 M h M c q 2 + 2m 2 M c (Mh + m 2 ^ 2 co 2 q 2 + Jjn 2 M c ^M h + m 2 )co 2 - m 2 M^M c S2 2 
+ m 2 M c (M h + m 2 )L 2 A 2 (L 2 )« 2 - m 2 (k h + m 2 )A 2 (L 2 )(m 1 L 1 - m^O 2 ^ 

- 2m 1 m 2 M c ? 1 cOjq 1 - Im-^nigMeCo 2 + m^rngMcLj A^Lj^fi 2 

- A l( L l)( m l L l - m 2 L 2)^ 2 ]qi + m 2 M c C (yh + fix h) + m 2 M c K y h = 0 (A54) 

By using the results of equation (A49) , equation (A43) can be written as the x h 
equation: 

M c M g (x h - 2f2y h - « 2 x^ + Mg^mjLj - m 2 L 2 )S2 2 - 2m 1 M s flq 1 + 2m 2 M s fiq 2 


+ ^M s + M c )C(x h - Oy h J + ^M s + M c )Kx h = 0 


(A55) 
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It should be noted that this equation contains a steady internally applied force propor- 
tional to £2 2 which disappears when the centrifuge is balanced, that is, nijLj = m 2 L 2 
and by using equations (A50), (A53), and (A54), equation (A45) can be expressed as 
follows: 

M c M h M s (y h + 2£2x h - £2 2 y h ) - 2m 1 M c M s C 1 w 1 q 1 - jm^McMgCO 2 + m 1 M c M s L 1 A-j^L j^£2 2 
- mjMg A^LjJ^mjLi - mgLj^^jq^ + 2m 2 M c M s £ 2 a; 2 q 2 + |” 1 2 M c M s a) 2 
+ m 2 M c M s L 2 A 2 (L 2 )£2 2 - m 2 M s A 2 (L 2 )(m 1 L 1 - m 2 L 2 )s/jq 2 + M c (m s + M h )C(y h + £2x h ) 

+ M C (M S + M h )Ky h = 0 (A56) 

Equations (A49), (A50), (A53), (A54), (A55), and (A56) represent the governing equations 
of motion for the configuration and are summarized in the following section subject to the 
restriction that m^Lj - m 2 L 2 = 0. 


Summary of Reduced Equations 

It should be recalled that in developing the equations of motion, the centrifuge was 
constrained to rotate at a constant angular velocity £2 relative to the space station. The 
space station was assumed to be constrained against rotational motion. In addition, non- 
linear terms were neglected. In summarizing the final dimensional form of the equations 
of motion, one additional assumption is made, that is, it is assumed that the centrifuge in 
the undeformed state is statically balanced about the hub as defined by the following 
equation: 

m^Lj - m 2 L 2 = 0 (A57) 


The equations of motion may now be summarized as follows: 

*s - 2S2y s - S2 2 x s - j^(i h - «y h ) - £ x h = 0 


(A58a) 


y s + 2SJx s - n 2 y s - ^(y h + S2x h ) - ^ y h = 0 


(A58b) 


9 Mo + M r /. N M s + M c 2mj£2 , 2m 2 £2 . 

- 2£2y h - £2 2 x h + „ c(x h - £2y J + Kx h - q x + q 2 = 0 


h M S M C 


M S M C 


M r 


(A58c) 
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.. . 9 1V1 «! lvl h /• \ M + M, 2m 1 

y h ♦ 2S* h - n 2y h ♦ C(y h + fi* h ) + Ky h L±A ,, 

m l"l m l L l a l( L l) f H 2m 2 ? 2 w 2- I“2"2 m 2 L 2 A 2( L 2) n 1 

_ M h + M h J q l + M h q 2 + _ M h + M h j q 2 " 


ixijWj j£2 


qi 


( M h + m l 

H 

M h 



2m ? C 2 a)„ m„w 2 m 9 L 9 A„fL 9 )o 2 l r . K v 

- -TT^ + ■ 2 2 J V 2; q 9 - — fy h + S2xJ) - JL y = o (A58e) 

M h 2 |_ M h M h J* 2 M h V h M h v JOfc!; 


( M h + m 2 y 

)4 

M h 



M h 




- O 2 + 


( M h + m 2) L 2 A 2( L 2) n 


M h 1 M 


q l + * nX h) + y h = 0 (A58f) 
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CENTRIFUGAL STIFFENING EFFECT 

In this appendix the equation of motion for the forced response of a rotating canti- 
lever beam is derived. The lateral deflection is approximated by a single term in a 
modal-form expansion. The resulting equation of motion is shown to contain a term pro- 
portional to the square of the rotational rate which includes the centrifugal stiffening 
effect. The tip of the cantilever beam is shown to follow a parabolic path. 

Symbols 

a(x,t) axial shortening, equation (Bl) 

EI(x) bending stiffness 

f(x,t) distributive loading 

L length of beam 

M generalized mass, equations (B4) 

m concentrated tip mass 

m(x) distributive mass 

Q(t) generalized force, equations (B4) 

q( ) generalized coordinate, equation (B2) 

t time 

X,Y rotating reference frame 

x coordinates along X-axis 

u(x,t) lateral deflection parallel to Y-axis 

X dummy variable for x 
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A(x) defined in equations (B4) 

6(x - L) Dirac delta function at x = L 
c p(x ) assumed mode shape 

£2 rotational rate of beam 

a) natural frequency of nonrotating beam (eqs. (B4)) 

u> n natural frequency of rotating beam (eq. (B14)) 

A dot over a variable indicates differentiation with respect to time. 

A prime over a variable indicates differentiation with respect to x. 

Analysis of Centrifugal Stiffening 

Consider the rotating centilever beam as shown in sketch (a) where £2 is the 
angular rate of rotation. 


u(x,t) 

- X 

Sketch (a) 

For any station x and time t there is an axial shortening a(x,t) due to a lateral 
deflection u(x,t). Thus, a point at coordinate x in the undeformed state is shifted 
toward the origin by the amount a(x,t) where 



a(x,t) = i C u ,2 (X,t) dX 
l ^0 


(Bl) 


If Lagrange's equation for this system is used and a single term in a modal form 
expansion for the deflection u(x,t) is assumed, that is 


u(x,t) = 0(x) q(t) 


(B2) 


the resulting equation of motion, the higher order nonlinear terms being neglected, is as 
follows: 


Mq + 


m(o^ - £2^) + £2^ f xm(x) A(x) dx 
J 0 


q = Q(t) 


(B3) 
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where 


~\ 


M = f $^(x) m(x) dx 
J 0 

Q(t) = f f(x,t) 0(x) dx 
J 0 

" 2 = sio EI & , " (x 3 2,te 

A(x)=J [0’(X)] 2 dX 


(B4) 




Equation (B3) is a forced undamped form of equation (A58e) of appendix A for a fixed hub 
support. Note that by the use of equation (B2) and the last relation of equations (B4), 
equation (Bl) can be written as 


a(x,t) = | A(x) q 2 (t) 

For the beam with a concentrated tip mass, that is, 

m(x) = m5(x - L) 

where 5(x - L) is the Dirac delta function, equation (B3) becomes 

q + |a> 2 + ft 2 [L A(L) - ljj q = Q(t) 


where 


A(L)=J L [0’(X)] 2 dA 


(B5) 


(B6) 


(B7) 


(B8) 


By further stipulating a uniform cantilevered beam and using a mode shape 0(x) which 
corresponds to that for a concentrated static load at the tip, 

EI(x) = El 


0(x) = -~(3Lx 2 - x 3 ) 
2L 3 


(B9) 


from which is obtained 


M = m 


= v 
ilS > 


ml 


L A(L) = | 


(BIO) 


J 
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Substituting equations (BIO) into the equation of motion (B7) yields 

mq + |m(u>2 - + 6m j^2jq _ (Bll) 

where S2 2 j s the centrifugal stiffening term. Thus, 

5 

q + (« 2 +|n 2 )q = l Q (t) (B12) 

where Q(t) is the generalized force associated with the mode 0(x) as defined in equa- 
tions (B4). Equation (B12) is a second-order linear differential equation with constant 
coefficients. The general solution will depend on the form of Q(t); however, the homoge- 
neous form of the equation is 

q + (to 2 + | 0 2 )q = 0 (B13) 

which is oscillatory and has a natural frequency o> n 

co n = co 2 + | O 2 (B14) 

1 9 

Thus, the dependence of u> n on the rotational rate O by virtue of the — £1* term 

5 

includes a centrifugal stiffening effect. (See refs. 1, 2, and 3.) 

For a beam with a concentrated tip mass, the desired form of equation (B5) becomes 

a(L,t) = I A(L) q 2 (t) (B15) 

which describes a parabola with axis on the undeformed beam axis, regardless of the 
assumed mode shape 0(x). However, the coefficient of the parabola A(L)/2 is related 
to the assumed mode shape and determines the amount of shortening and hence stiffening 
that may result. An aspect showing the importance of this shortening or stiffening may 
be seen in equation (B3) where the motion is unstable when the stiffening is insufficient, 
that is, 

S2 2 f x m(x) A(x) dx < Mf2 2 
J 0 

For the mode shape assumed in equation (B9), the stiffness is adequate and equation (B15) 
becomes 

a(L,t) = — q 2 (t) (B16) 

5L 

Because the tip motion for the cantilevered beam traverses a parabola which may be 
represented by a circle for small deflections, an alternate configuration as shown in 
appendix C was developed. 



APPENDIX C 


EFFECTIVE HINGE POINT 

In appendix B the equation of motion for a rotating cantilever beam was deter- 
mined in order to assess the centrifugal stiffening effect. In this appendix an alternative 
rotating system as shown in sketch (b) is analyzed. When this system is made equivalent 
to the rotating cantilever beam of appendix B, the effective hinge point may be easily 
defined. 


Symbols 

torsion spring simulating beam bending 
L undeformed length of two-link mechanism 

m tip mass 

q lateral displacement coordinate 

X,Y rotating coordinate axes 

x position along X-axis 

a hinge-point location relative to length L 

A composite function of length and hinge-point location (see eq. (C5)) 

6 coordinate of rotation about hinge point 

£2 rotational speed 

o o uncoupled natural frequency about hinge point (see eq. (C4)) 

Dots over symbols denote time derivatives. 

Analysis of Rotating System 

Consider the two-link rotating system shown in sketch (b). The system is con- 
strained to rotate uniformly at a constant rate £2. Coupling of the two links is provided 
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by means of a torsional spring of stiffness lq. At any time, the deflection of the mass m 
from the X-axis may be determined in terms of 6 . The net movement of m, it should 
be noted, consists of a motion transverse to the X-axis and a motion along the X-axis. If 
the equation of motion for the system is derived and then linearized by using a small angle 
assumption on 6 (sin 6 ~ 6 ), the resulting equation is: 


m(L - aL )^6 + maLO^(L - aL)0 + kj.0 = 0 
If the lateral deflection is defined as 

q = (L - aL )6 

then equation (Cl) may be written as 

■>2 


L - aL, 


q = 0 


(Cl) 

(C2) 


(C3) 


where the substitution 


a, 2 = 




m(L - aL) 2 

has been made. 

Comparing equations (C3) and (B7) for unforced motion and defining A = A(L) 

gives 


(C4) 


or 


A = 


a = 


1 

L - aL 

(C5) 

LA - 1 
LA 

(C6) 


where 0 g a S 1 and aL locates an effective hinge point as measured from the actual 
center of rotation (x = 0). Thus, equation (C6) defines the normalized position of a cen- 
ter about which the tip of a beam effectively pivots during small deflections. 
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Imaginary, j co* 



Figure 4.- Locus of roots as function of rotational speed for even stability polynomial in p-plane and p 2 -plane. Rotational speeds Q* increase from 0 for a through g. 



Normalized natural frequency. 



Normalized rotational speed. 


a* 


Figure 5.- Uncoupled loci of stability determinant. 
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Normalized rotational speed, SI* 


Natural frequencies of basic model in rotating reference frame. Mi = M 2 = ^3 = °* 01 » = 

aj = Q 2 = 1/6. Asymptotes shown in phantom. 



o* 

o 



Gondola masses Gondola masses 

deflect equally do not deflect 

q | = q 2 = Constant q| r q 2 = 0 


(a) Lower critical-speed mode. 


(b) Upper critical-speed mode. 


Figure 7.- Upper and lower critical-speed modes, m^ = m 2 ~ m. 








4,= -100 0 = 0 ° 





“>*, = 0 W n| = 0 

(Unstable) (Stable) 



* * * ^ 

w n,= 0 OJn } = 0 

(Unstable) (Stable) 



<t)* 2 = I 000 5=0° 






w * 3 = 1.822 0 = 39 . 5 ° 




a » nz =.828 0 = 108 . 7 ° 




Figure 10.- Typical mode shapes for symmetric undamped model with oscillatory instability. 
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< u * ( = 2.258 0 = 127 . 5 ° 
(Unstable) 


& = 3.8 



w * 2 = 2.258 0 = 127 . 5 ° 



**,= 2.766 0 = 123 . 6 ° 
(Unstable) 



**, = 3.105 0 = 121 . 7 ° 
(Unstable) 


< i ?„ 2 = 2.766 0 = 123 . 6 ° 




a .„ 2 = 3.105 0 = 121 . 7 ° 



Figure 10.- Concluded. 






"*3 = 6.003 0 = 75 . 0 ° 
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Uncoupled mode 



_L 

3 


4 


5 


Normalized rotational speed, Cl* 

Figure 11,- Natural frequencies in rotating frame with pj = p 2 = °.3, M 3 = 0.01, u;* = = 0.1 and aj = a 2 = 1/6. 

Asymptotes shown in phantom. 



&=0 
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aVi 
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(Unstable) (stable) 



* * - 

w n| = 0 td n | = 0 

(Unstable) (Stable) 






w* 2 = (-192 0 =45 3° 

w n 2 = J. 277 0 =49.3° 




w* 2 = 1.585 0 = 56.8° 



a, * 2 = 1.814 0 = 59.5° 



a;*, = 4.01 1 0 = 112.2° 


0 ? 


, = 5.976 0 = 75.3° 



0 = 90.0° 



0 = 90.0° 


wn 3 =l.58l 0=0° 



0,^=1592 0 = 5 . 6 ° 



w* n3 =1.802 0 = 25.0° 


w* 3 = 1.880 0 = 28 7° 



"* 3 = '964 0 = 32.1° 


w* 3 = 2.236 0 = 40.3° 



w *n 3 = 2.429 0=44.5° 



oj* 3 = 7.1 91 0 = 62.6° 



Figure 12.- Typical mode shapes for symmetric undamped model without oscillatory instability. 
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Normalized rotational speed, Cl* 

Figure 14.- Result of inboard effective hinge location on natural frequencies in rotating reference frame. 
Mx - M 2 = 0 . 1 , M 3 = 0 . 01 , and cdf = v*= 0 . 1 . 










Normalized rotationc 



I I I ! ! ! 1 

0 .2 .4 .6 .8 1.0 1.2 


Stiffness ratio, k ( /K 


Figure 20.- Effect of heavy centrifuge arm damping on the oscillatory divergent instability regions. M 3 = 0.01, cq = a 2 = l/ 6 , 

q = Z > 2 = °* 5 - and C = °‘ 01 - 
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